FRACTIONAL ^-DIFFERENCE EQUATIONS 
ARISING FROM THE CALCULUS OF VARIATIONS 



RUI A. C. FERREIRA AND DELFIM F. M. TORRES 

Abstract. The recent theory of fractional ^-difference equations introduced 
in [N. R. O. Bastos, R. A. C. Ferreira, D. F. M. Torres: Discrete-time 
fractional variational problems, Signal Process. 91 (2011), no. 3, 513-524], 
is enriched with useful tools for the explicit solution of discrete equations 
involving left and right fractional difference operators. New results for the 
right fractional h sum are proved. Illustrative examples show the effectiveness 
of the obtained results in solving fractional discrete Euler-Lagrange equations. 



1. Introduction 

The fractional calculus is a generalization of (integer order) differential calcu- 
lus, allowing to define derivatives (and integrals) of real or complex order [251 130) . 
It is a mathematical subject that has proved to be very useful in applied fields 
such as economics, engineering, and physics (3J [2lJ [22] [29] . Several definitions of 
fractional derivatives, including Riemann-Liouville, Caputo, Riesz, Riesz-Caputo, 
Weyl, Grunwald-Letnikov, Hadamard, and Chen derivatives, are available in the lit- 
erature [21[T31[I7J[53]. The most common used fractional derivative is the Riemann- 
Liouville [TJ [13j [TBI HZ] • Analogously, one can define a discrete fractional derivative 
in different ways. In 1989, Miller and Ross introduced the discrete analogue of 
the Riemann-Liouville fractional derivative and proved some properties of the frac- 
tional difference operator [24]. More results on the theory introduced by Miller 
and Ross are given in the works of Atici and Eloe IS]- See also [TUJ [TTJ [T5] . 
and [7] for applications to the Gompertz fractional difference equation and tumor 
growth models. Regarding other fractional discrete definitions, we refer the reader 
to [U HU [28] and references therein. Here we follow [9], i.e., we adopt a more 
general fractional ^-difference Riemann-Liouville operator. The presence of the 
h parameter is particularly interesting from the numerical point of view, because 
when h tends to zero the solutions of the fractional difference equations can be seen 
as approximations to the solutions of corresponding Riemann-Liouville fractional 
differential equations 0113] (cf. Proposition [3]). 
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In the recent work of Bastos et al. [9], necessary optimality conditions of first 
and second order are proved for the fractional /i-diffcrence variational problem 

(1) C ^ = E L Vl?h(th)), ( a A a h y) (th), ( h A?y) (th)) h — > min , 

(y(a) = A), (y(b) = B), 

as well as transversality conditions when the boundary conditions y(a) = A or 
y(b) = B are not given (see Section [5] for definitions and notations). The main 
result of [9] gives an Euler-Lagrange type equation for problem ([1]), but no clues 
are devised for the solution of such fractional /i-difference equations. Instead, some 
examples are solved numerically |9] . Here we develop further the subject of the 
calculus of variations within the fractional discrete setting, by obtaining explicit 
solutions to the fractional difference Euler-Lagrange equations [HI [H] . 

Our results are given in Section [2] where we prove some new formulas for the 
fractional ^-difference operator. The obtained results are then used in Section [3] to 
solve two illustrative examples of (fTJ) , for which the global minimizers are explicitly 
found in exact form. This is in contrast with [5J [S], where all the solutions are 
obtained via approximated numerical computations. 

2. Main Results 

Before stating and proving our results, we introduce some definitions and nota- 
tions. Let h > and put T = {a, a + h,a + 2h, . . . , b} with net and b = a + kh 
for k £ {2,3,...}. Let us denote by Tf the set of real valued functions defined on 
T, <Jh(t) =t + h, and Ph(t) = t — h. 

Definition 1. For a function / s J~t the forward /i-difference operator is defined 
as 

(A h f)(t) = f^h(t))-f(t) tG { aa + ha + 2h p ( b) y 
h 

while the /i-difference sum is given by 

(.AfcV)(*) = J2 f ( kh 1 h ' t£{a,a + h,a + 2h,..., a h (b)}. 

*=* 

Definition 2. For arbitrary x,y £ M. the /i-factorial function is defined by 

x (v) _ h v r(f + i) 
h ~ m + l-y)' 

where T is the Euler gamma function. We use the convention that division at a 
pole yields zero. 

In [5] it is remarked that in the case h = 1, then coincides with the falling 

(v) 

factorial power. One also expects to see that converges to x v when h tends to 
zero. Since this is not addressed in [S], we prove it here. 

Proposition 3. For x > and y £ K ; 

(2) lim x { ? ] = x v . 

/i-X) 
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Proof. Equality (j2J is a straightforward consequence of the following well-known 
asymptotic formula for the Gamma function: 

(see, e.g., inequality (33) and Corollary 3 in [23]). Indeed, starting from the defini- 
tion of a:^ and introducing the new variable t = x/h + 1 — y, we have 

.(») _ x y nt + y) 



X 



ll 



(t + y- i)y r(t) 



for any y G K. We obtain ([2]) taking the limit t — > +oo or, equivalently, the limit 
h -> 0+.' □ 

The motivation for the next definition can be found in [9]. 

Definition 4 ([9]). Let / G J^. The Ze/i and rzg/ii fractional h-sum of order 
1/ > are, respectively, the operators A ft " : .Ft — >• Fj.+ and TiA^" : Ft — >• Fj-, 

= {a ± a ± i'/i ± h, ■ ■ ■ , b ± vK\, given by 



- — V 



( q a-<7) (*) = £ (* - ^(^))r i} /(^)^ 



^=7 

I 

77 



( A A 6 ""/) (t) = ^ £ (kh-a h (t)) ( r 1} f(kh)h. 



We define ( a A°/) (t) = /(t) and ( h Ag/) (t) = /(*). 

Definition 5. Let < a < 1 and set 7 = 1 — a. The left fractional h-difference 
Q A^/ and the right fractional h-difference hA^ f of order a of a function / G Ft 
are defined, respectively, by 

(„A£/) (t) = (A/jqA^ 7 /) (t) , tefa + TM + THft,..., p h (6) + 7 M , 

( h Af /) (t) = -(A hh Aff(t) , i G {a - 7 A, a - 7/1 - ft, . . . , p h (b) - 1 h}. 

Remark 6. We define fractional sums/differences for functions on a bounded do- 
main. This is done so, because of the problems of the calculus of variations we con- 
sider here. Nevertheless, one can use our definitions for functions with unbounded 
domains: an unbounded domain from above for the left fractional sum/difference, 
an unbounded domain from below for the right fractional sum/difference. 

Let us now recall a result that will be used later in finding solutions to the 
boundary value problems originated from the fractional /i-difference calculus of 
variations. 

Theorem 7 (Theorem 2.10 of [9]). Let f G .Ft and v > 0. Then, 

( a A^A h f) (t) = (A ha A^f) (t) (t - a) { r 1] !{a) 

for all t G {a + uh, a + vh + h, . . . , phip) + vti). 

The next lemma permits to shorten the proofs of our main results. Essentially, 
it allow us to borrow information from the formulas obtained in [5]. 
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Lemma 8. Let f £ Jy and v > 0. Then, 

( a A f ;7)(i)=/ i "( f A^m/(s/ i )) , tefo + ^a + ^ + ^.J + Wi}. 
Proof. We have 

( Q A^/) (t) = ^ ^ (t - ^(fc/O^/W 



f 7 — y 



fe=4 



and therefore the proof is done. □ 

Lemma 9 (Lemma 2.3 of [5]). Let fj,,v be two real numbers such that fj,,fi + v £ 
E\{...,-2,-l}. Then, 

(„Ar, -> w = r^+t+V ^ ' ' = M + M + " + h • • • 

Corollary 10. Suppose that £, £ + i/ e . . , -2, -1}. TTien, 

( A- S ^ ( S - a + (t) = r^|^i) (* - a + 

/or a?/ i £ {a + a + vh + h, . . .}. 

Proof. The result is a simple consequence of Lemma [5] and Lemma [§] □ 
Corollary 11. Suppose that + v £ R\{. . . , -2, -1}. Then, 

( h A^s ^ (6 - M - .)<"*>) ( t ) = pgjffif^ (& — M — t)'"*^ 

/or all t £ {b — vh, b — vh — h, . . .}. 
Proof. Using Corollary [101 we get 

_& 

^A 6 -" a .->(6-M-a)^ ) )(*) = I 4r £ (fcfe - ^(t))^ (6 -/i - fe/i)^ h 
J (kh + b-n-a h (t))^(-kht^h 



I» 



b — f_i-t 



-pA,^ -> 4^) {b — fi — t) = r( _g* ^ (& - M - C 
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and the proof is complete. □ 
We now state and prove the law of exponents for the fractional /i-difference sums. 
Theorem 12. Let f G and [i,v > 0. Then, 

(3) ( a+vh A^ a A^f) (t) = ( a A-^f) (t), 
where t G {a + (/i + u)h, a + (/i + v)h + h, ■ ■ ■ , b + (fi + v)h}; and 

(4) ( h A-^ hh A^f) (t) = ( h A- ( ^f) (t), 
where t G {b — (// + 6 — (/i + — /i, . . . , a — (fi + 

Proof. We prove ([3]) only, ((4]) being accomplished analogously. First, note that if 
u = 0ovfi = u = 0, then the equality is valid by definition. Therefore, assume that 
f-l,y-l+/i€ K\{. . . , -2, -1}. Then, 



a+vh A^ a A^f){t) 

e (*-^(^))i^ x) e^-^m)^/^) 



i»r( M ) 
r(z/)r(u) 



,2 T!— (/*+") ^- r-l-M 

r^fXM) E E (t-afcK. + r + lJfc])^^)^/^) 

r ^ s v 1 



h 



-77 



E ((„-!)* V* -> 1} ) (f - a*(rfc))/(rfc) 



-(/-H-f) 



" E r^(*-^))r /, " 1) /(rfc)=(.A^)/)(t), 



i» r(^ + M ) 

which shows the intended equality. □ 

The next theorem is crucial in order to solve some fractional ^.-difference Euler- 
Lagrange equations (see the examples in Section [3]). 

Theorem 13. Let < a < 1 and f G J-j. Then, 

(5) ( h Ag7) (t) = 0, te{a-(l-a)/i,a-(l-a)/i + /i,...,p fc (6)-(l-a)/i}, 
i/ and cm/?/ i/ 

(6) /(t) = ^( & _(i_ a )/ l _f)^- 1 ) ) te{a,a + fc,...,6}, 

r(a) 

where c is an arbitrary constant. 
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Proof. If / is given as in (j6]), we immediately get ([5]) using Corollary [TT] Suppose 
now that equality in ([5]) holds in the mentioned domain. Then, by definition of 
fractional difference, 
(7) 

f/,A^ (1_a ° /)(*)= c, t e {a-(l-a)h,a + h-(l-a)h,...,b-(l-a)h}, cel. 



Applying the operator h^ b ^ 1 _ a ^ h to both sides of equality in (ff|. and using (j4]) of 

Theorem HI we get (^A" 1 /) (<) - c ( h A-_? (1 _ a)fc l) (*), te{fl-fc,o ) ... ) 6- ft}. 
Corollary [Til now implies that 

(8) (hAt 1 /) (*) = =7 C (b-(l-a)ft * 6 {a -ft, a,..., b -ft}. 
1 (a + 1) 

An application of the operator A^ to both sides of the equality in (J5J) gives 



/K(*)) - c - (6 - (1 - a)h - <J h (t))[ \ t€{a-h,a,...,b- 2ft}, 
1 (a + 1) 



or 



(9) /( t )=*(&_(i_ a )fc ie{a,a + ft,...,6-ft}. 



Setting t = b — h in © we get f(b)h — ch a , i.e., equality in © is also valid when 
t = b. □ 



Remark 14. Similar steps as those done in the proof of Theorem 1131 permit us to 
prove the following equivalence: for < a < 1, c G R, and / G J-t, we have 



(foA£7) (i) = c, i e {a - (1 - a)ft, a + ft- - (1 - a)ft, . . . , b - h - (1 - a)h] 
if and only if 



„, . c(b + ah — ba — a 2 h — t) + da , . , ,( a -i) r , , ^ 

f(t) = - r(a+l) (1 1)}" 1J , te{a,a + ft,...,fe}, 
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where d is an arbitrary constant. Indeed, from Corollary ITTI and Theorem W2\ we 
have 



ft 



A£7)(t) = c^ ( h A b - {1 - a ^f)(t) = -ct + d 



( h A~ (1 ~ a) /) (t)=c(b + ah-t-b- ah) + d 
( h A~ (1 ~ a) /) (*) = c(b + ah-t)- c(b + ah) + d 
(^Mi-a).^^ 1 "'/) (*) = c (,A-« (1 _ Q)/iS ^ 6 + aft - S ) (t) 
-[ c (6 + ^)-d](,A-« (1 _ Q);i l) (i) 

- U V/) (*) - + ^ - - c( r{ a tV d (6 " (1 " a)h ~ t)[a) 

^ f(t) = 2 r(^f^ + l^_ _ c(b + ah)-d _ (a _ x) 

p / b-(l- a )/ t -t + A 

«*•/(*) = ?^ — r(& + a/i-t)/i a - x ' '' ' ' 



T(a + iy ' ^' r ^-(i-#-t + 1 _ (a _ 1) j 

c(6 + ah-ba- a 2 h - t) + da 
= r(a + 1) (& — (1 — o;)/i — t)X ; - 

We end this section enunciating the analogue of Theoremll3lfor the left fractional 
/i-difference. 

Theorem 15. Let < a < 1 and f G Jy. Then, 

(«A£/ ) (*) = 0, i G {a + (1 - Qi)/i, a + (1 - a)h + h, . . . , p h (b) + (1 - a)h}, 
if and only if 

/(*) = T ^(t-(l-a)h- a) ( «- 1] , t G {a, a + h, . . . , b}, 
where c is an arbitrary constant. 

Proof. The proof is analogous to the one of Theorem [T3] □ 

3. Applications to the Calculus of Variations 

We now give two examples of application of our results. The main achievement 
is to obtain explicit solutions for some problems of the calculus of variations. In 
this section we omit the subscript h in ah and p^- For convenience of notation we 
write y a {t) =y{a(t)). 

Example 16. Let us consider the following data: let a G R, h > 0, b = a + kh with 
k G {2,3,.. .}, and < a < 1. Moreover, let A and B be to given real numbers. 
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We want to find a function y G Tt that solves the problem 

(10) C(y)= J2 UA%y) 2 (th)h ^ min, y(a)=A, y(6)=fl. 

By [SI Theorem 3.5] we have that if y is a minimizer of C given in (|10p. then 

(11) ( fc A£ (6)o Ag0) (t) = 0, t6{o,o + fc,...,i-2h}. 

Remark 17. At a first glance the sum in (|10[) and the equation in (|lip seem to 
be meaningless due to the possible values of the variable t. However, they aren't 
by the fact that the authors in |9] used the following notation for the difference 
operators: 

( a A%f) (t) = ( a A a h f) (t + (1 - a)h), te{a,a + h,...,b-h}, 
( h Af /) (t) - ( h A?/) (t - (1 - a)A), i G {a, a + h, . . . ,b — h}. 

An application of our Theorem [13] to the equality in (jlip gives 

( Q A£y) (t) = ^ (p(6) - (1 - a)h t)^ , t e {a, a + h, . . . , p(6)}, 
with c € K or 

(A^A^y) (t) = ^ (p(b) - (t + (1 - a)ft))^- 1} , t € {a, o + fe, . . . , 

We now remember Remark [T71 and apply the operator a +(i-a)hAj^ a to both sides 
of this equality. From Theorems [7] and [T2] it follows that 

(12) 

( a+i ^ a)h A^A ha Al- a y) (i) = ^ ( a+(1 _ a)h A^« S i ^ (p(6) - s)^) (t) 

° y(^) = f^y G+u-^Ar* ^ M 6 ) - 4 a_1) ) (*) 

+ ^(*-(a + (l-a)/i))l a - 1) y(a), 
r(a) 

with £ G {a + /i, a + 2/i, . . . , 6}. The constant c is determined by the end condition 
y(b) = B. 

Remark 18. We point out that if a = 1 we get the "straight line" connecting the 
points (a, A) and (b, B) as the solution of the Euler-Lagrange equation (TT21 . i.e., 
y(t) = fff (i - a) + A This result can be found, e.g., in [121 120]. 

We now show that the function y given by (j!2l) furnishes in fact a global min- 
imum to the problem (1101) . To do that, we recall the fractional /i-summation by 
parts formula obtained by the authors in [S] (we continue to use here the notation 
mentioned in Remark 1 17)) . 

Theorem 19 (Theorem 3.2 of [9]). Let f and g be real valued functions defined on 
{a, a + h, . . . , b — h} and {a, a + h, . . . ,6} , respectively. Fix < a < 1 and put 
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7 = 1 — a. Then, 

E (° A ^) (th)h = hy(p(b))g(b)-hy(a)g(a)+Y^ (h*% b) f) {th)g°{th)h 

*=T7 *=T 

ft, h. 

73(a) ^^j, , „a7-i) f/ , U j, V" M, i „^^(7-i). 



T(7+l) 



E(i« + 7^-«),V E (^ + 7^-^(«))i 7_1) /(^)ft 



Before proceeding, we need the following definition: 

Definition 20. We say that a Lagrangian L(t, u, v) : T x R 2 — > R is jointly convex 
in (u, v) if 

L(i, w, u) — L(t, u , v) > (u — u')L u (t, u , v ) + (v — v')L v (t, u , v), u, u , v, v G R, 

provided the partial derivatives L u and L„ exist. 

We are now able to prove the following theorem. 

Theorem 21. Consider the set S = {/ G J-~t : /(a) = A,f(b) — B}. Suppose that 
the Lagrangian L(t, u, v) : T x R 2 — > R o/ i/ie minimization problem 

is jointly convex in (u,v). Assume that the function y G S satisfies the Euler- 
Lagrange equation for this problem, i.e., 

(13) L u [y](t) + (h^ a p(b) L v [y]) (t) = 0, te{a,a + h,...,b-2h}, 

where [y](s) = (s,y l7 (s), ( A^y) (s))). Then, y furnishes a global minimum to C in 
the set S . 

Proof. Let y G S be an arbitrary function. Suppose that y G S satisfies equation in 
(1131) . Since L(t, u, v) is jointly convex in (u, v), we get, with the use of Theorem ll91 

E {L[y\ - Mm h>J2 i(y a - v a )Lu[y] + { a Kv - a Ky)L v [y]}h 

= E fa - y a ) ( L «[y] + ^ {b) L v m) } h 



The theorem is proved. □ 

It is clear that the Lagrangian L(t,u,v) = v 2 in (|10|l is jointly convex in (u, v). 
Therefore, the function y defined in (|12j) furnishes a global minimum to ([T0|) . 



We end solving another fractional difference problem of the calculus of variations. 
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Example 22. Let a g R, h > 0, b = a + kh with k € {2, 3, . . .}, < a < 1, and A 
and B be two given real numbers. We consider the following variational problem: 



(14) £(!/) = E 



-i ^ 

±( a AZy) 2 (th)-y°(th) 



y(a)=A, y(b) = B. 



The Euler-Lagrange equation for problem (fT4")l is 

(15) ( A A£ (6) „Agy) (f) = 1, i 6 {a, a + /!,.--,&- 2/i} . 

In view of Remark [T4j we get from equality in fj 15[) that 

, .„ . ,. p(b) + ah — p(b)a — a 2 h — t + da , ,,. , . , N f Q -ii 

( a A&) (*) = — Yd + 1) ~ (1 ~ a)h ~ t)h ' 

for a constant d £ K to be determined. Following the same steps as those done for 
Example UM we get 

Ma\ rM ( a —a , b- p{b)a-a 2 h~ s + da \(a-i)\ M 

(16) = I a+ (i_ a ) h A h rTaTTj yp(°)- a )h ) (*) 

+ f ^y(i-(«+(l-")^), ( r 1) 2/(«). 

for f G {a+/i, a+2/i, . . . , fe}. Finally, we show that the Lagrangian L(i, u, w) = v 2 — u 
is jointly convex in (u, v). Indeed, for u, w, m', h'e! we have 

v 2 - u - (i/ 2 - u') > -(u - v!) + (v - v')2v' &{v- v 1 ) 2 > 0. 

We conclude that y given by (TT6|) is the global minimizcr of (|14p . 
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